A randomized polynomial-time algorithm 
for the Spanning Hypertree Problem on 3-uniform hypergraphs 
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Consider the problem of determining whether there exists a spanning hypertree in a given fc-uniform 
hypergraph. This problem is trivially in P for k = 2, and is NP-complete for fc > 4, whereas for 
fc = 3, there exists a polynomial-time algorithm based on Lovasz' theory of polymatroid matching. 

Here we give a completely different, randomized polynomial-time algorithm in the case A; = 3. The 
main ingredients are a PfafRan formula by Vaintrob and one of the authors (G.M.) for a polynomial 
that enumerates spanning hypertrees with some signs, and a lemma on the number of roots of 
polynomials over a finite field. 



INTRODUCTION 

A (finite) hypergraph G = (V, E) consists of a finite set 
V (the vertex set) and a set E of subsets of V (the hy- 
peredges), each of cardinality at least 2. We also write 
V{G) = V and E{G) = E. When all the hyperedges have 
the same cardinality k, we say that the hypergraph is k- 
uniform. In the case fc = 2 we are dealing with ordinary 
(simple) graphs. We say that H C G is a, sub-hypergraph 
if V{H) C V{G) and E{H) C E{G); a sub-hypergraph is 
spanning if V{H) — V{G). A hypergraph H is a hyper- 
tree if it is connected and there are no cyclic sequences 
of vertices and hyperedges 



Vi,Ai,V2,A2, 



,vi,Ae,vi+i (ve+i = vi) 



such that £ > 2, ut ^ Vi-^i and Vi, Vi+i G At. An example 
of a 3-uniform hypergraph, together with a spanning sub- 
hypergraph that is a hypertree, is shown in Figure [TJ 

We will deal here with issues in Computational Com- 
plexity Theory [5|. An introduction to the subject which 



Figure 1: A 3-uniform hypergraph G, and a spanning hyper- 
tree T C G, whose edges are in darker gray. 




includes the class RP of probabilistic polynomial-time 
problems, pertinent to this paper, can be found in Chap- 
ters 2-4 of Talbot and Welsh [15|. 

In this paper we are concerned with the following de- 
cision problem: 

fc-UNiFORM Spanning Hypertree (fc-SHT): 
Given a fc-uniform hypergraph G — (V, E), deter- 
mine whether there exists a spanning hypertree. 

Of course, every connected graph contains a spanning 
tree, so fc-SHT is trivially in P for fc = 2 (it suffices to 
check whether G is connected). On the other hand, for 
fc > 3 it is not true that every connected fc-uniform hy- 
pergraph contains a spanning hypertree, and the decision 
problem is highly nontrivial. 

Our main result here is to provide an RP algorithm for 
the fc-Uniform Spanning Hypertree problem when fc — 
3. After a first version of this paper was completed, we 
have learnt from Andras Sebo that there is actually a 
polynomial-time algorithm for this problem, coming as a 
specialization of Lovasz' algorithm for matching on linear 
2-polymatroids [1, Q. However, Lovasz' techniques are 
completely different from ours, and we believe that our 
more algebraic approach is of independent interest. We 
remark that for fc > 4 the spanning hypertree problem is 
NP-complete by a result of C. Thomassen which appears 
in [3, Theorem 4]. (We thank Marc Noy for bringing 
this argument to our attention.) Moreover, the same 
argument shows that the corresponding counting problem 
is ttP-complete already for fc > 3. We will briefly review 
Thomassen's argument at the end of this introduction. 

Organization of the paper. The bulk of this paper has 
two parts. In the first part, we discuss the main ingredi- 
ent of our RP algorithm, namely the Pfaffian-Tree Theo- 
rem of [m which expresses a signed version of the multi- 
variate spanning-tree generating function of a 3-uniform 
hypergraph as a Pfaffian. Then in the second part, we 
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describe our algorithm, first intuitively, and then more 
formally, and sketch the analysis of time- and space- 
complexity. This part is fairly standard in complexity 
theory, but we hope that the partly expository presen- 
tation of the various concepts involved will be useful for 
the interdisciplinary audience (such as ourselves) we have 
in mind. Finally, we end the paper with some specula- 
tions and directions for further research suggested by our 
work. 

To conclude this introduction, here is, then, Thom- 
assen's argument showing that k-SHT is NP-complctc for 
> 4. 

We recall that an exact cover in a hypergraph G = 
(y, E) is a subset E' <Z E oi the hyperedges such that 
every vertex of G belongs to exactly one hyperedge in 
E' . (In the special case where G is an ordinary graph, 
an exact cover is nothing but a perfect matching.) Now 
consider the following decision problem: 

Exact cover by /c-sets (XfcC): 
Given a A:-uniform hypergraph G = (V, E), deter- 
mine whether there exists an exact cover. 

X3C is known to be NP-complete, and is classified as 
problem |SP2] in Garey and Johnson (It is NP- 

complete even when restricted to 3-partite hypergraphs, 
in this case being called 3-Dimensional Matching (3DM, 
|SP1]).) Implicitly, analogous statements hold as well 
for any fc > 3. Conversely, X2C is polynomial, by 
matching techniques, even in its optimization variant, 
e.g. by Gallai- Edmonds algorithm (see [l^). On the 
other hand, the corresponding counting problem (equiva- 
lently: counting perfect matchings on arbitrary graphs), 
is known to be {jP-complete (VaHant [16]). 

Now, given an arbitrary fc-uniform hypergraph G, 
Thomassen constructs a (fc + l)-uniform hypergraph G' 
as follows: add an extra vertex and let 

E{G') = {e\J {i<}\e ^ E{G)} . 

The key observation is that spanning hypertrees of G' 
correspond bijectively to exact covers of G (with the ob- 
vious bijection, namely deleting * from each hyperedge). 
Thus, any algorithm for (fc -1- 1)-SHT provides an al- 
gorithm for XfcC. In other words, XfcC is reducible to 
(fc + 1)-SHT. 

From what we said above about XfcC, it follows that 
fc-SHT is NP-complete for fc > 4, and counting spanning 
hypertrees in a 3-uniform hypergraph is jJP-complete, as 
asserted. 



A PFAFFIAN FORMULA 



on G is 



Zg{w) 



E 



n 



WA 



(1) 



TeT{G) AeE{T) 



where T(G') is the set of spanning hypertrees of G, and 
the WA are commuting indcterminates, one for each hy- 
peredge A G E{G). 

Assume that G is fc-uniform. If G has a spanning hy- 
pertree, then necessarily 



(2) 



where n is the number of hyperedges in each spanning 
hypertree. Therefore we will assume ([2]) from now on. 

It will be convenient to consider G as a sub-hypergraph 
of 1C{N, fc), the complete fc-uniform hypergraph on N ver- 
tices, which has an hyperedge for all unordered fc-uples 
A ^ C [N] := {1, 2, . . . , iV}. We denote 

Zic{N,k)iw) by Zn^k{w), and T{IC{N,k)) by Tn,k- Note 
that the degree of Zn^kiw) is n, as given in ([2|). The 
spanning hypertree generating function Zq{w) of an ar- 
bitrary fc-uniform hypergraph G can be obtained from 
Zn,k{w) by setting the weights wa to zero for hyperedges 
A not in G. 

A classical result by Kirchhoff is that, for k — 2, the ex- 
pression Zn.2{w) (and therefore Zg{w) for any graph G) 
is given by a determinant. Defining the N x N Laplacian 
matrix L as 



-w,j i ^ j 



(3) 



and taking a whatever {N — l)-dimensional principal mi- 
nor L{i(j) (i.e. with row and column removed - and 
remark that iV— l = niffc = 2), one has 

Theorem 1 (Matrix- Tree) 

^n,2(w) =detL(io). (4) 

As is well-known, Kirchhoff's formula allows one to 
count spanning trees on a graph G: putting wa = 1 if 
A G E{G) and wa = otherwise, the determinant @ 
gives the cardinality of the set T(G) of spanning trees 
on G. For later use, we remark that this formula also 
shows that counting spanning trees on graphs is in P, as 
determinants can be evaluated in polynomial time. 

For general fc-uniform hypergraphs, no such formula is 
known if fc > 3. Moreover, a determinantal expression 
for Zn^kiw) is unlikely to exist, since counting spanning 
hypertrees is ftP-complete if fc > 3, as discussed in the 
introduction. 

We remark in passing that, nevertheless, the number of 
spanning hypertrees on the complete hypergraph 1C{N, fc) 
is known: generalizing the classical result by Cayley for 
fc = 2, one has that 



Let G = (y, E) be a finite hypergraph on N vertices. The 
multivariate generating function for spanning hypertrees 
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as can be found in Q (see also references therein). 

From now on, we will consider the case fc = 3. Om' 
starting point is a recent result, due to A. Vaintrob and 
one of the authors (G.M.) ll|, which states that an alter- 
nating sign version of the spanning hypertree generating 
function Zn^3{w) is given by a Pfafhan. 

We will denote this modified polynomial by Z*^{w). 
For a given ordering of the vertices of the hypergraph (or, 
equivalently, a labeling of the vertices with integers from 
1 to N), a sign function e(T) for hypertrees T is defined 
(more details are given later); then 



; AeE{T) 



(6) 



To state the result, let e^fe be the totally antisymmetric 
tensor (i.e. Cijk = if two or more indices are equal, 
^ijk — 1 a i < j < k OT any other cyclic permutation, 
and eijk = — 1 if i < A; < j or any other cyclic permuta- 
tion) and define a A'^-dimensional antisymmetric matrix 
A, with off-diagonal elements 



Ay - ^ ^ijkW^i^j^k} ■ 

Then one has, for any principal minor A(io): 
Theorem 2 (Pfaffian-Hypertree 



(7) 



(8) 



(The original formulation in uses indeterminates yijk 
which are antisymmetric in their indices. The corre- 
spondence with our notation here is simply that yijk — 

In order to give meaning to equations ([6]) and ([8]), we 
must define the sign function e(T) : 7^^3 — > {±1}. Sev- 
eral equivalent definitions exist, and all of them require 
making some arbitrary choices; the proof that the result- 
ing sign is actually independent of these choices (devel- 
oped to full extent in 11|) can be performed inductively 
in tree size, studying the invariances in the elementary 
step of adding a 'leaf edge to a tree. 

It is worth stressing, however, that the precise deter- 
mination of this sign function is not actually used in 
our algorithm (and not even in our proofs of complex- 
ity bounds). 



The first definition of t{T) given in is as follows. 
For S C [N] , call ts the permutation which rotates cycli- 
cally the elements of S (in their natural order) , and keeps 
fixed the others. Then, for a given ordered n-uple of hy- 
peredges (Ai, . . . , A.^) forming a tree T, define the per- 
mutation f = ■ ■ ■ TA„ ■ This permutation is composed 
of a single cycle of length N. It is thus conjugated to the 
"canonical" A^-cycle T[jv] , i.e. there exists a such that 



Then actually the signature e{a) does not depend on the 
ordering of the hyperedges, but only on T, and taking 
e(r) :— e{a) is a valid definition, and the appropriate 
one for ([8]) to hold. 

An equivalent definition is as follows. Consider a pla- 
nar embedding of the tree, in such a way that for each hy- 
peredge Aa, if {ia,ja,ka) denote its three vertices cycli- 
cally ordered in the clockwise order given by the em- 
bedding, one has eijk = 1- (Such an embedding always 
exists.) Then construct the string p of 3n symbols in [N] 
(recall that n is the number of hyperedges) correspond- 
ing to the sequence of vertices visited by a clockwise path 
surrounding the tree, starting from an arbitrary vertex. 
Each of the N vertices occurs in this string, but some ver- 
tices appear more than once. Now remove entries from 
this string until each of the N vertices appears exactly 
once, thus getting a string p oi N distinct elements of 
[N]. If we interpret this string as a permutation, then 
e{T) — e{p), despite of the arbitrariness of the choices for 
the planar embedding, the starting point and the extra 
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a T\ 



\N] ' 



(9) 



entries we choose to remove. An example is given in 
Figure 3.1]. 

Instead of using a planar embedding, one can also 
describe this procedure by choosing a given vertex as 
root (say, io), and orienting the edges accordingly (so 
that each edge has a "tip" and two "tail" vertices, and 
all edges are oriented towards the root). Now the ver- 
tices inside any hyperedge Aa can be uniquely ordered 
{ia,ja, ka) such that ia is the tip and ei^j^k^ = 1- With 
this notation, the sign e{T) coincides with the sign of a 
certain monomial in the iV-dimensional exterior algebra 
(or real Grassmann Algebra): 

e(io)A(/\e(i„)Ae(A:,)) = e(r) e(l) A • ■ • Ae(iV) . (10) 



In order to see the equivalence with the previous defi- 
nition of e(T), it suffices to take an appropriate planar 
embedding of T and to observe that we can perform the 
removal of extra entries in such a way that in the final 
string p, the two tails of each hyperedge always come 
consecutively. The ordering of the hyperedges in the 
product on the L.H.S. of (fTO|l is irrelevant at sight, as 
an hyperedge with odd size has an even number of tails, 
which thus are commuting expressions in exterior alge- 
bra. Nonetheless the true invariance is stronger, as it 
concerns also the choice of the root vertex, and much 
more, as implied by the preceeding paragraphs. The ex- 
pression pop above is an efficient way of computing e{T) 
for a given tree T. 

This last definition for e(T) has also the advantage of 
driving us easily towards a proof using Grassmann vari- 
ables of formula ([8]). In fact, for both equations ([4]) and 
([5]), say with iq = 1, we can recognize the expressions for 
Gaussian integrals of Grassmann variables, "complex" 
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and "real" in the two cases respectively 

Zr,^2{w) = I V{^,i:) ^,^,e^^'^^^^-^^^^^'-^^^ ■ (11) 
Z* 3(w) = / V{e) 6*1 e5:'-'..'^'(«'«j+«^«'=+»'^'»0 . (12) 



It is combinatorially clear that such an integral will gen- 
erate an expansion in terms of spanning subgraphs, with 
one component being a rooted tree and all the others 
being unicyclic. Then one realizes that, because of the 
anticommutation of the variables, unicyclics get contri- 
butions with opposite signs, which ultimately cancel out. 
However this mechanisrn^occurs in slightly different ways 
in the two cases (cfr. [4] and respectively for more 
details). 

Remark. The original proof of the Pfaffian-hypertree 
theorem in 11| was by induction on the number of hy- 
peredges using a contraction-deletion formula. Another 
proof was given in [l3| which exploits the knot-theoretical 
context which originally lead to the discovery of the for- 
mula. The proof using Grassmann variables alluded to 
above is due to Abdesselam [1] . Finally, yet another proof 
was given by Hirschman and Reiner 6] using the concept 
of a sign-reversing involution. 



TOWARDS A RANDOMIZED 
POLYNOMIAL-TIME ALGORITHM 

Equation ([6]) provides us with a multivariate alternating- 
sign generating function, which is a polynomial, and iden- 
tically vanishing if and only if our hypergraph G has no 
spanning hypertrees. Can we use this to decide efficiently 
whether G has a spanning hypertree? 

Recall that in the classical Matrix- Tree theorem for or- 
dinary graphs there are no signs, and upon putting edge 
weights equal to 1 for edges in the graph, and equal to 
zero otherwise, we get the number of spanning trees on 
our graph. Moreover, calculating numerically the deter- 
minant of a matrix can be performed by Gaussian elimi- 
nation in polynomial time. 

Let us apply the same idea to the Pfaffian-Hypertree 
formula Since det A(io) = (Pf A(io))^, we can again 
evaluate Z* ^{w) in polynomial time, if, as before, we set 
the hyperedge weights equal to 1 for hyperedges in a 3- 
uniform hypergraph G, and equal to zero otherwise. But 
because of the signs in Z* ^, this is not the number of 
spanning hypertrees in G! 

One way out is to evaluate det A(io) at some random 
set of numerical weights (but keeping weights equal to 
zero for hyperedges not in the hypergraph). If one gets a 
non-zero result, this would certainly prove that the mul- 
tivariate generating function is non-zero, and hence pro- 
vide a certificate of the fact that a spanning hypertree 
exists. Conversely, if many evaluations at random in- 
dependent points are zero, one starts believing that the 



graph has no hypertrees at all. This naive idea can be 
formalized within the framework of the RP complexity 
class. 

In complexity theory, the class of Randomized Poly- 
nomial-time problems (RP) contains problems for which, 
given any instance, a polynomial-time probabilistic algo- 
rithm can be called an arbitrary number of times in such 
a way that: 

• If the correct answer is 'False', it always returns 
'False'; 

• If the correct answer is 'True', then it returns 'False' 
for the t-th query with a probability at most 1/2, 
regardless of the previous query results. 

We now discuss in more detail how to construct an RP- 
algorithm from the Pfaffian-Hypertree formula ([8|) . Note 
that, of course, one has P C RP C NP. 



GAUSSIAN ALGORITHM OVER FINITE FIELDS 

Calculating numerically the determinant of a matrix 
using Gauss elimination is commonly thought to be of 
polynomial time complexity (at most cubic). If this is 
certainly true for "float" numbers (but suffers from nu- 
merical approximations), some remark is in order for "ex- 
act" calculations. Indeed, in this case we have to choose 
a field, such as Q, which is suitable for exact numerical 
computation through a sequence of sums, products and 
inverses, which, in the complexity estimate above, have 
been considered an "unity of complexity" (a variant is 
possible, in which one works in the ring Z and recursively 
factors out a number of g.c.d.'s). This is however not an 
innocent assumption. For example, if one works with ra- 
tional numbers with both numerators and denominators 
having a bounded number of digits d, in general, during 
the Gauss procedure, one may suffer from an exponen- 
tial growth of this number (as the l.c.m. of two d-digit 
integers may well have 2d digits). 

An improved choice, and which makes the analysis sim- 
pler, is to consider finite fields GFg, for which the com- 
plexity of operations is uniformly bounded. A primer in 
finite fields can be found, for example, in the textbook 0] . 
The field GF, exists for q a power of a prime, p'^, as a 
quotient of a set of polynomials with coefficients in Zp 
by a polynomial of degree h, irreducible in Zp. Also for 
finite fields, polynomials of degree n have at most n dis- 
tinct roots. Recall that = iff a = in any field, so 
that also in GF^ one has det A = only if Pf A = 0. 

Two specially used cases are q prime, in which GFq 
just coincides with Z^, equipped with the product mod- 
ulo q, and q a power of 2, which is mostly used in Coding 
Theory. In both cases, the arithmetic operations -|-(-, •), 
x(-,-) and l/(-) are performed in a very efficient way. 
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if a 'Space' requirement of order q is allowed for stor- 
ing a table of discrete logarithms (this will turn out to 
be subleading for our problem). As a result, if a given 
value for q is chosen, calculating the determinant of a 
TV-dimensional matrix in GF((7) takes a time of order 
N^^logq)"^ (improved to N^logq if a fast preprocessing, 
concerning the table of logarithms, is performed 0, US])- 
Remark that for us it is important to control the depen- 
dence not just on N but also on q, because in our context 
q will depend on N. 



ROOTS OF POLYNOMIALS OVER A FINITE 
FIELD 

The analysis of the previous sections naturally induces us 
to study the roots of polynomials over Galois fields GF^. 
In particular we need an upper bound as follows: 

Lemma 1 Let f{xi,...,Xn) be a non-zero polynomial 
in n variables of total degree d with coefficients in the 
finite field GFq. Then f has at most dq^~^ roots. 
In other words, the probability that a randomly chosen 
(ai, . . . ,a„) G GF^ is a root of f is < d/q. 

For the convenience of the reader, we include a proof 
of this lemma. An essentially equivalent proof formu- 
lated in the language of probability theory is given in 
[l5l Thm. 4.2]. 

Recall that the field GF^ has q elements. We say that 
(ai, . . . , On) G GF^ is a non-root of / if /(ai, . . . , a„) ^ 0. 
We must show that / has at least {q — d)q'^^^ non-roots. 

The proof is by induction on n, the number of vari- 
ables. For n = 1, a non-zero one- variable polynomial of 
degree d has at most d roots. Now assume n>2. Write 

m 

f{xi,...,Xn)=^x\gi{x2,.-.,Xn) , 

i=Q 

where gm and m < d. Observe that gm has degree 
at most d — m. By the induction hypothesis applied to 
gm, we know that gm has at most {d — m)q"'~'^ roots, and 
therefore we can find at least {q — d + m)q"~'^ non-roots 
of gm- For every (02, . . . , a„) which is a non-root of gm, 
consider 



THE ALGORITHM 

The lemma [1] stated above implies that, for any 3-uniform 
hypergraph G with 2n + I vertices, if one evaluates 
(Z* 3)^ at some values {waIas-ECG) random uniformly 
sampled from GF^ (and, of course, wa = if A ^ E{G)), 
one gets always zero if G has no spanning hypertrees, and 
obtains zero although G has some spanning hypertrees, 
with probability at most n/q. If q > 2n, we are within 
the framework of the RP complexity class. 

So our algorithm is just as follows. Given G, choose 
a root vertex and an ordering (once and for ever), and 
build the corresponding matrix A of indeterminates as in 
([7]) (as a matrix of lists of edge-labels, with signs). Build 
the table of GF^-logarithms for an appropriate value of 
q. Then, for a given fault tolerance e = 2^*^, repeat k 
times the following probabilistic algorithm: 

• extract the values wa independently identically dis- 
tributed in GFq, 

• evaluate A in numerical form for these values; 

• evaluate the determinant numerically by Gauss 
elimination (in GF^); 

• if the result is non-zero, return "there are trees" , 
and break the cycle. 

Then, if the cycle terminates without breaking, return 
"there are no trees with probability 1 — 2^^" . 

If we consider a generic value of g > n, when the cycle 
terminates without breaking we know that there are no 
trees with probability at least 1 — (n/q)''. For a given 
value of e, we thus get an upper bound on the time com- 
plexity (up to a multiplicative constant) 



If both n,e^^ — > 00, optimization in q suggests to take 
q ^ n/e and perform a single query, (instead of taking 
q — 2n and performing 0{ — loge) queries). This saves an 
extra factor min(logrt, — loge), and gives a complexity of 
order max(log7i, — loge). 

PERSPECTIVES 



f{xi,a2, . . . ,a„) . 

This is a non-trivial polynomial in xi of degree exactly m, 
therefore it has at most m roots and at least q — m non- 
roots. Thus we have found at least {q — d-\-m)q'^~^{q — m) 
non-roots of / (namely {q — d -\- m)q^~'^ possibilities for 
(02, . . . , a„) times q~ni possibilities for ai.) It is easy to 
check that {q- d + m)q''-^{q ~ m) > {q - d)q"-^. This 
completes the proof. □ 



Given the datum of a 3-uniform hypergraph G = (V, E), 
a positive integer s, and a set of integer- valued edge costs 
{cyi}, consider the cost function for spanning hypertrees 
'H{T) — J2AeE{T) Then, one has the decision prob- 
lem of determining if G has any spanning hypertree T 
with 7i{T) < s. This problem is essentially equivalent to 
the corresponding optimization problem, of finding the 
spanning hypertree of minimum cost, and has a complex- 
ity at least as large as the one of S-SHT. At the time of 
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this writing, we don't know whether Lovasz' polymatroid 
matching techniques can answer this problem, which is 
not addressed in [1, [§1 . Be that as it may, it would also 
be interesting to understand if an RP algorithm exists 
for this problem, through an extension of the technique 
described in the present paper. 

Another issue worth investigating may be how to use 
the exact expression for the sign e{T) and maximize the 
fraction of terms getting the same sign in ((6]), in order to 
use the Pfaffian-Hypertree formula in all of its strength. 
This may lead to a notion of Pfafhan orientation for 3- 
uniform hypergraphs, analogous to the notion of Pfaffian 
orientation (or Kasteleyn orientation) for graphs [Ts] . In- 
deed, our story is somewhat parallel to the story for per- 
fect matchings (and also Lovasz algorithm is a general- 
ization of the Edmonds-Gallai algorithm for matchings 
on graphs). In particular, if we construct a 3- uniform 
hypergraph G" from a graph G using the procedure de- 
scribed in the introduction, the Ffafhan-Hypertree for- 
mula ([8]) for G" becomes the well-known alternating sign 
multivariate generating function for perfect matchings of 
G given by a Pfaffian. Now in the graph case, there is the 
notion of Kasteleyn orientation which serves to get rid of 
the signs in the Pfaffian and allows therefore to count 
perfect matchings in this way. A Kasteleyn orientation 
exists for planar graphs, but not in general. Is there an 
interesting class of "Pfaffian orientable" 3-uniform hyper- 
graphs for which the Pfaffian-Hypertree formula can be 
used to count spanning hypertrees exactly? 
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